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Abstract 

We propose a one-dimensional rice-pile model which connects the ID BTW 

sandpile model ( Phys. Rev. A 38, 364 (1988) ) and the Oslo rice-pile model 

( Phys. Rev. lett. 77, 107 (1997) ) in a continuous manner. We found that 

for a sufficiently large system, there is a sharp transition between the trivial 

critical behaviour of the ID BTW model and the self-organized critical (SOC) 

behaviour. When there is SOC , the model belongs to a known universality 

class with the avalanche exponent r = 1.53. 

PACS numbers: 64.60.Lx. 05.40-hj, 64.60.Ht. 05.70.Ln 
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Bak et al proposed "self-organized criticality" (SOC) as a framework to understand the 
dynamics of driven, dissipative systems These systems, via self-organization, reach the 
steady state which is characterized by a power-law distribution of the sizes of avalanches. 
The authors of Ref. originally used a cellular automata, now referred to as the BTW 
sandpile model, to illustrate their ideas. Moreover the SOC behavior was also found in 
some biological |^ and economical [§] systems. Different variations of the BTW sandpile 
were proposed and studied The behaviors of these sandpiles are somewhat different 

depending on whether the rules of evolution are based on the absolute sand heights of the 
pile the local slopes @J^, or the Laplacians of the sand height function 0. Besides 

the numerical simulations many different methods were also used to treat the SOC prob- 
lems. Dynamical mean field theory [H gives a unified description of some stochastic SOC 



systems including the BTW sandpile model and the forest fire model [10|. Langevin-type 



approaches [|ri[] have been used on a phenomenological basis. Furthermore, a real space 
renormalization group method [|T^ provided good estimates of the avalanche exponents. Fi- 
nally, it has been shown |^ that a large class of sandpiles were Abelian and this property 
leads to a particularly simple equiprobable partitioning in configuration space that allows 



to extract some exact results. Numerical simulations of high dimensional BTW model 
were recently performed to determine the upper critical dimension where the avalanche dis- 
tributions are characterized by the mean-field exponents. The idea of SOC also stimulated 



much interest in the granular matter and some experiments [|I4| were done to investigated 
whether real sandpile display SOC behavior. In order to make comparison with theories and 
models, a group of researchers in Oslo did experiments on real rice piles [0 and showed 
that under some conditions a real rice pile displays SOC behaviors. In fact, for grains with 
a large aspect ratio the system self-organizes into a critical state. They explained this result 
with the increased friction and different packing possibilities. By measuring the transit time 
after the pile has reached the stationary critical state, they found that the distribution of 
the transit times follows the form: 
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P{T,L) = L-f^F{T/U), (1) 

where T is the transit time, L is the system size. The scahng function F{x) is constant for 
small X and decays as power law with a slope a = 2.4 ± 0.2 for larger x. 

To take into account the changes in the local slopes observed in the rice pile experiment, 
Christensen et al proposed a rice-pile model, hereafter called Oslo model [p!6HT8[| , where the 
critical slope for each site is a dynamical variable. The Oslo model is based on a linear array 
of cells labelled by i, where i = 1,2, and an integer variable h{i) assigned to each of 
them, with a wall at i = and an open boundary at i = L + 1. Here h{i) is called the 
local height of the rice pile at site i. The local slope assigned to each site i is defined as 
z{i) = h{i) — h{i + 1) for i = 1,2, L. Initially, the system is empty, i.e., h{i) = \fi. The 
dynamics of the model consists of deposition and relaxation, and the relaxation process is 
considered to be fast compared to the deposition time scale. At each time step one grain 
is added to a column i which increase the height of i by 1, i.e., h{i) h(i) + 1 . With 
the dropping of rice grains, a rice pile is built up. Whenever there is active column, i.e., 
z{i) > z'^{i), where z'^{i) is a slope threshold, one grain of rice will be transferred from this 
column to its right neighbor, h{i) — > h{i) — 1 and h{i + l) h{i + l) + 1, and all the unstable 
sites topple in parallel. The critical slope z"^ of a site remains unchanged if the site is stable 
but assumes a new value 1 or 2 randomly every time a rice grain on this site has toppled. 
The toppling of one or more sites is called an avalanche event, and during the avalanche no 
grains are added to the pile. The avalanche stops when the system reaches a stable state 
with z{i) < z'^{i) Vz. The internal randomness in the critical slopes makes the Oslo model 
different from the ID BTW model. For the Oslo model, with arbitrary initial conditions, 
the system reaches a stationary state where the avalanche sizes are power-law distributed, 
while for the ID BTW model the size distribution of avalanche is generally not power-law. 
Since the only difference between the ID BTW model and the Oslo model is the presence of 
internal disorder in the latter, it is natural to consider that the criticality in the Oslo model 
is induced by the disorder ( randomness in the critical slopes) . One motivation of this letter 
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is to investigate the transition between the ID BTW model and the Oslo model. 

We modify the Oslo model as follows: When z{i) < z'^{i) Vz, the pile is stable and there 
isn't diffusion of particles. If, at site i, z{i) > z^{i), where z^{i) takes randomly a value 
1 or 2, then this site topples with a probability which depends on its slope z{i), namely: 
if 1 < z{i) < 2, the site i topples with a probability pi and if z{i) > 2, it topples with 
probability p2. Notice that if we set pi = and p2 = 1, the model becomes the BTW model 
with the critical slope Zc = 2, and if we set pi = p2 = 1, it is just the Oslo model. In a 
previous paper [jl9l, one of us generalized the Oslo model in a different way, in which the 



critical local slope can assume r different values, from 1 to r, where r is an integer. It is 
possible to show that the exponents for the avalanche-size and transit-time distributions 
are insensitive to the level of medium disorder (different values of r > 1) when the grain is 
dropped at a fixed position. So in the present model we will only consider the case r = 2, 
i.e., the critical slope takes randomly the value 1 or 2. In this study, we will restrict ourselves 
to the case where a grain is added to the site i=l. When a grain is dropped on the left-end 
site of the pile, it may make the site unstable. The site topples and transfers a grain of rice 
to its right neighbour and so on. And in this way avalanches occur. As in the literature, 
we define the size of an avalanche as the number of toppling. In this letter, we will take p2 
equal to 1 and pi < p2 because higher is the slope, higher is the jump probability. So by 
varying pi from to 1, we can change the model from the ID BTW sandpile model to the 
Oslo model in a continuous manner. 

We have performed extensive numerical simulations and investigated the effect of pi on 
the behaviors of the model. 

Let us first study the transport properties of the model. The transit time of a grain is 
defined as the time it spent in the pile, and the time is measured in the unit of additions of 
grains. When a grain slips out of the pile, we can measure its transit time T = Tout — Tin-, 
where Tout and Tj„ denote the input and output time of the grain. For the case pi = (the 
ID BTW model) when the stationary state is reached every newly-added grain will slip out 
of the pile instantly, thus the transit time is T = 0, and the average transportation velocity 



of grains, defined as < v >= L/ < T >, is infinite for this case. On tlie otlier liand, for the 
case pi = 1 (the Oslo model) previous studies show that the average transportation velocity 
scales with the system size as < f >oc L''' with 7 = 0.30 ± 0.10, indicating that < v >^ 
in the limit of infinite system size. It would be interesting to see how < v > changes from 
00 to when pi is varied from to 1. In fact, we found that there is sharp transition at 
Pi = 0. The general behavior of < > as a function of pi is shown in Fig. ||. One can see 
from Fig. Q that for a given system size L, the velocity < v > becomes constant when pi is 
larger than some value, say pf, where pi itself is dependent on the system size and becomes 
closer to when L becomes larger. The numerical results make us to consider that pi —>■ 
as L ^ 00. . When pi 0, < v >= 1, independent of the system size. So there is a sharp 
transition from < v >= 00 for pi = to < >= 1 for pi = 0^. 

At first think, the sharp transition may be surprising. But it can be understood by the 
following argument. It is clear that when pi is exactly 0, no newly-added grain will stay 
in the pile as long the stationary state is reached. So T = for every grain, and hence 
< T >= 0. For pi = 0+ however, the situation is different, and there is the possibility 
that some grains will be buried in the surface layer of the pile. These grains will stay in 
the pile for a very long time. Once they slip out of the pile, these grains, although very 
few in number, will make a significant contribution to < T > since their transit times are 
extremely large. It is the existence of these grains that makes < T > assume finite value for 
Pi = 0^. Thus the sharp transition here is induced by the tiny disorder. Between pi = O"*" 
and pi = p\i there is crossover behavior of < >, which is due to finite size effects. Since 
we expect — when L 00, we can also expect that for infinite system the transition 
takes place at pi = form < v >= 00 to < f >= 0. 

In what follows we shall show that the critical behaviors of the model belong to the same 
universality class of the Oslo model When pi is larger than pi. In Fig. 0, we plot the average 
velocity as a function of the system size for different values of pi greater than pi. It is clear 
that for large enough system {L > 100), the average velocity < v > scales as L^"^ , where 
7 ^ 0.23. Therefore the average velocity decreases with the system size, which is due to the 
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increase in the active zone depth with system size, as explained by Christensen et al . 



We have also studied the avalanche size and the transit time distributions for different 



values of the probability pi > p^. In a previous studies on the rice pile models |[T^ it was 
shown that the avalanche size distribution is of the form: 

P{S,L) = S-^G{S/L^), (2) 

with r = 1.53 ± 0.05 and D = 2.20 ± 0.05. We report in Fig. §.a our simulation data 
for the distribution of avalanche size for different values of system size. The distribution 
is a power law with the presence of a peak close to the cutoff size Sc oc L^. This peak is 
due to the finite-size effects which leads the system into a supercritical state, followed by a 
massive avalanche. Since on average each site must topple exactly once during an avalanche 
to transport a grain out of the pile when the stationary state is reached, the average size 



of avalanche is thus < S >= L. And this leads to the scaling relation D{2 — t) = 1 |T8 
The numerical results of the exponents, namely r ?a 1.53 and D ^ 2.20 are in quite good 
agreement with the scaling relation. The transit-time distribution can be described by the 
scaling form (1) with the same exponents as the Oslo Model {pi = 1), namely, f3 ~ 1.25 and 
u 1.25. The power-law exponent a for the large transit time is obtained as a ~ 2.4. ( 
See Fig. ^ for an example ). Fig. |^.a shows the avalanche-size distribution for a pile of size 
L = 400 and for several values of jumping probability pi. It is clear that the size exponent 
T is insensitive to the value of the probability pi greater than a critical value pi- Fig. ^.b 
gives the corresponding transit time distribution, which is nearly a constant for small transit 
time, and decays as a power law for larger transit time. As in Fig. |[a, the exponents remain 
the same for several values of pi greater than p^. Therefore, we can conclude that for pi 
greater than pi the size and transit time exponents are the same as the Oslo rice-pile model, 
thus the SOC state in our model is not affected by the fact that some grains topple with 
different jumping probability pi, when pi > pi- 

The crossover behaviors were also investigated. We collect statistics on the size of the 
avalanche and transit time of the grains for a value of pi very close to zero {pi = 10~^). As 
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illustrated in Fig. ^a, the avalanche- size distribution does not exhibit a power-law behavior. 
The data can be described by the following form: 

P(5,L) =L-^G(5/L") (3) 

which is confirmed by the good data collapse obtained with the exponents x = 0.94 and 
a = 0.94 and where the scaling function G{x) decays exponentially ( Fig. ^.b ). The 
form of Eq.(|^) is very different from the distribution form for the case pi = 0, which is 
P{S,L) = 6{S — L), where 6{x) is the Dirac Function. This confirms our statement that 
there is a sharp transition at pi = 0. In Fig. |^.a, the transit time distribution is plotted. It 
is visually apparent that the distribution peaks at the vicinity of the system size. A data 
collapse is shown in Fig. |^.b, which is obtained by rescaling the original data according to 
Eq. (|l|) with f3 = 1.2 and u = 1.2. The exponents used here are a little smaller than that 
for pi > p\. We expect that when pi ^ the exponents v and j3 shall become even smaller 
and finally approach 1. In this case the scaling for average velocity < v >~ L'"* = L^~^ will 
give < f >~ L*^ = 1, which is in consistent with the results that < >= 1 for finite system 
at pi = 0'^. 

In summary, we have investigated a one- dimensional model of a rice pile where the sites 
with higher slopes have more chance to topple (with a probability p2 = 1) while the sites 
with lower slopes topple with a probability pi < P2- Depending on the value of pi, our 
model exhibits three behaviors: the trivial critical behavior for pi = (ID BTW model), 
the crossover behavior due to the finite size of the rice-pile and the self-organized critical 
behaviour for pi greater than a certain critical value pi. In fact, for a sufficiently large 
pile, pI goes to zero and therefore, there is a sharp transition between the trivial behaviour 
and the SOC behaviour at pi = 0. When the system exhibits the SOC behaviour, the 
exponents did not depend on the value of the probability pi and our model belongs to the 
same universality class as the Oslo rice-pile model. 
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FIG. 1. The average velocity < v >, defined in the text, as a function of the probabihty pi 
for several values of the system size. 
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Fig. 2 




FIG. 2. The average velocity < v > as a function of the system size L, for two values of the 
probability pi. 
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Fig. 3 
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FIG. 3. (a) Log-Log plot of the avalanche size distributions for several values of the system 
size L with pi = 0.8. (b) Data collapse of the curves displayed in (a) according to Eq. (|2[) with 
the exponents r = 1.53, D = 2.2. 
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Fig. 4 




FIG. 4. (a) Log-Log plot of the transit-time distribution for several values of the system size 
with pi = 0.8. (b) Data collapse of the curves displayed in (a) according to Eq. (||). The best fit 
to the numerical data gives the slope a = 2.4 
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Fig. 5 
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FIG. 5. (a) Log-Log plot of the avalanche size distribution. The best fit gives the slope 
r = 1.53 (b) Log-Log plot of the transit time distribution. The best fit gives the slope a = 2.4 In 
both (a) and (b) the system size is L = 400. 
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Fig. 6 
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FIG. 6. (a) Log-Log plot of the avalanche size distribution with pi = 0.01. (b) Semi-log plot of 
the date rescaled according to Eq. (^. The curve is quite straight line, indicating the exponential 
form of the scaling function G{x). 
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Fig. 7 
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FIG. 7. (a) Log-Log plot of the transit-time distribution with pi = 0.01. (b) Log-Log plot of 

the rescaled data with the exponents /3 = 1.20 and ly = 1.20. 
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